Coherent QCD multiple scattering in proton-nucleus collisions 
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We argue that high energy proton-nucleus (p+A) collisions provide an excellent laboratory for 
studying nuclear size enhanced parton multiple scattering where power corrections to the leading 
twist perturbative QCD factorization approach can be systematically computed. We identify and 
resum these corrections and calculate the centrality and rapidity dependent nuclear suppression 
of single and double inclusive hadron production at moderate transverse momenta. We demon- 
strate that both spectra and dihadron correlations in p+A reactions are sensitive measures of such 
dynamical nuclear attenuation effects. 



PACS numbers: 12.38.Cy; 12.39.St; 24.85.+p 

Copious experimental data pj from central Au+Au re- 
actions at the Relativistic Heavy Ion Collider (RHIC) 
has generated tremendous excitement by pointing to the 
possible creation of a deconfined state of QCD with en- 
ergy density as high as 100 times normal nuclear mat- 
ter density @- In order to diagnose its properties 0, 
we first need to quantitatively understand the multiple 
scattering between a partonic probe and the partons of 
the medium in simpler strongly interacting systems, for 
example p+A. The importance of such theoretical inves- 
tigations has been recently stressed by the measured in- 
triguing deviation Q of the moderate-pT single inclusive 
hadron spectra from the independent multiple nucleon- 
nucleon scattering limit [j| in d+Au reactions at RHIC. 

In this Letter we present a systematic calculation of 
the coherent multiple parton scattering with several nu- 
cleons in p+A collisions, which has so far been consid- 
ered a challenge for the perturbative QCD factorization 
approach jfj. Coherent multiple scattering and the well- 
studied elastic multiple scattering co-exist in nuclear col- 
lisions and their relative role depends on the probes (or 
observables) and the underlying dynamics. Elastic scat- 
tering has played an important role in understanding the 
/cT-broadening and the Cronin effect Q. It dominates 
when the probe is so localized that the quantum cor- 
relation between the different scattering centers can be 
neglected. In terms of the factorization approach, con- 
tributions from coherent multiple scattering to a physi- 
cal cross section are power suppressed in comparison to 
the leading hard partonic processes Q. However, as we 
demonstrate below, when enhanced by the large nuclear 
size these may become important at small and moderate 
transverse momenta at RHIC. 

Hard scattering in nuclear collisions requires one large 
momentum transfer Q ~ xP 3> ^qcd with parton mo- 
mentum fraction x and beam momentum P. A simple ex- 



ample, shown in Fig.^a), is the lepton-nucleus deeply in- 
elastic scattering (DIS). The effective longitudinal inter- 
action length probed by the virtual photon of momentum 
g M is characterized by 1/xP. If the momentum fraction 
of an active initial-state parton x <C x c = l/2mjv?'o ~ 0.1 
with nucleon mass mjy and radius r , it could cover sev- 
eral Lorentz contracted nuclcons of longitudinal size ~ 
2ro(mN I P) in a large nucleus 8]. In the photon- nucleus 
frame, Fig. ^a), the scattered parton of momentum I 
will interact coherently with the partons from different 
nucleons at the same impact parameter || . The on-shell 
condition for t fixes the initial-state parton's momentum 
fraction to the Bjorkcn variable x — xb — Q 2 /{2P ■ q) 
but additional final state scattering forces the rescaling 
x = x B (l + (A 1 / 3 - 1)/Q 2 ) @. Here, A is the nuclear 
atomic weight and the scale of power corrections 
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FIG. 1: Coherent multiple scattering of the struck parton 
in deeply inelastic scattering (a) and in proton-nucleus colli- 
sions (b). A set of resummed two gluon exchange diagrams (c) 
in p+A reactions. 
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with the matrix element (p\F 2 \p) = | lim x ^o xG{x, Q 2 )- 
In Ref. Q we demonstrated that the hard parton inter- 
actions for any number of coherent multiple scattering in 
DIS are infrared safe, which explicitly verifies the factor- 
ization theorem. For £ 2 = 0.09 — 0.12 GeV 2 the calcu- 
lated reduction in the DIS structure functions, known as 
nuclear shadowing, is consistent with the xb-, Q 2 - and 
A-dependence of the data and the small Q 2 modification 
of the QCD sum rules . Coherence should play a simi- 
lar role in virtual quark and gluon exchanges in hadronic 
collisions, as shown in Fig. 1(b). 

Before we resum the power corrections in p+A reac- 
tions and derive the modification to the single and dou- 
ble inclusive hadron production cross sections we need 
to address the validity of factorization in a hadronic en- 
vironment. The colored soft gluon interactions between 
the incoming and/or outg oing hadrons can potentially 
ruin the factorization jlOl Hl| . Nevertheless, factoriza- 
tion in hadronic collisions was proved to be valid at the 
leading power Q as well as the lea ding power corrections 
(D(l/Q 2 ) It has been argued |l3j] that factorization 
can be extended to the calculation of (A 1 / 3 / Q 2 ) N -type 
nuclear size enhanced power corrections in p+A collisions 
at all powers of N. 



In DIS, as shown in Fig. 0a), only diagrams with mul- 
tiple final-state interactions to the scattered parton lead 
to medium size enhanced power corrections [j. On the 
other hand, in p+A collisions, shown in Fig. 0b), all dia- 
grams with either final-state and/or initial-state multiple 
interactions could in principle contribute. However, the 
coherent multiple scattering is a result of an extended 
probe size, which is determined by the momentum ex- 
change of the hard collisions. As shown in Fig. 0b), 
once we fix the momentum fractions x a and z\ , the effec- 
tive interaction region is determined by the momentum 
exchange g M = (x a P a — P c j ' z\Y . In the head-on frame 
of q — P b , the scattered parton of momentum I inter- 
acts coherently with partons from different nuclcons at 
the same impact parameter. Interactions that have taken 
place between the partons from the nucleus and the in- 
coming parton of momentum x a P a and/or the outgoing 
parton of momentum P c /z\ at a different impact param- 
eter are much less coherent and actually dominated by 
the independent elastic scattering 

We pursue the analogy with the DIS results of Ref. 
and express the lowest order single and double inclusive 
hadron production cross sections as follows: 
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dx b S(x b - x b ) F ab ^ cd (x b ) , (2) 



dx b S(x b - x b ) F ab ^ cd (x b ) , (3) 



where ^2 abcd runs over all parton flavors, 4>a/N{x a ), 
4> b /i\f(x b ) are the parton distribution functions 
(PDFs) [13 and D hl/c (ztLD h2/d (z 2 ) are the frag- 
mentation functions (FFs) [ly. The dependence on the 
small momentum fraction x b is isolated in the function 
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where the 2 — > 2 on shell squared partonic matrix el- 
ements \M ab ^, cd \ 2 are given in 01 ■ I n Eq. 101 the in- 
variants T = -p Tl ^/Se~ y \ U = -p Tl VSe Vl and the 
momentum fraction from the nucleus x b — —T/(z\S + 
U/x a ). In Eq. © Aip = ip 2 - <fi, z 2 = Zip T2 /p T i, 
and the solution for the momentum fractions carried by 
the partons from the nucleon and the nucleus 
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spectively. In both Eqs. J5J and © the possible if-factor 
for the lowest order formula and all renormalization and 
factorization scale dependences are suppressed. Ratios of 
cross sections in p+A and p+p reactions are insensitive 
to their choice. 

As in the leading v4 1 / 3 -enhanced contribution to 



the cross sections comes from the coherent scatterings 
with N nucleons shown in Figs. 0b) and0c). The two 
gluon exchange with an individual nucleon gives rise to 
the characteristic scale of high twist corrections £ 2 = 
0.09 - 0.12 GeV 2 , defined in Eq. 0. For fixed N, the 
final state interactions of parton d, illustrated in Fig 1(c), 
replace the integral over the small x b dependent function 
F ab ^cd(x b ) in Eqs. and @ by 
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Here, the hard scale t = q 2 = (x a P a — P c jz\) 2 and C d 
is a color factor depending on the flavor of parton d. 
C q (q) — 1 and C g — Ca/Cf = 9/4 for quark (antiquark) 
and gluon, respectively. By resumming all contributions 
in Fig. 1(c), N — 0, • ■ ■ , oo in Eq. 0, we derive the 
nuclear modified cross sections for single and double in- 
clusive hadron production in p+A collisions. These have 
the same form as Eqs. and with the following 
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FIG. 2: Top: suppression of the single and double inclusive 
hadron production rates in d+Au reactions versus pt for ra- 
pidities j/i = 0,1,2,3 and 4. £ 2 = 0.12 GeV 2 . Data on for- 
ward rapidity h~ attenuation is from BRAHMS y|. Bottom: 
impact parameter dependence of the calculated nuclear mod- 
ification for central 6 = 3 fm, minimum bias &min.bias = 5.6 fm 
and peripheral b = 6.9 fm collisions. The trigger hadron 
Pti = 1.5 GeV, y\ = 3 and the associated hadron y 2 = 0. 
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Similarly to the DIS case, resumming the coherent scat- 
tering with multiple nucleons is equivalent to a shift of 
the momentum fraction of the active parton from the nu- 
cleus and leads to a net suppression of the cross sections. 
The i-dependence of this shift in Eq. © indicates that 
the attenuation increases in the forward rapidity region. 

We consider d+Au reactions at RHIC at y/S = 
200 GeV with the deuteron and the nucleus moving with 
rapidities y max = ±5.4, respectively. Small pr ± = 1 GeV 
and forward y± = 4 hadron production corresponds to nu- 
clear Xb > 1.3 x 10~ 4 where coherence and our calculated 
power corrections will certainly become important. Dy- 
namical nuclear effects in multi-particle production can 
be studied through the ratio 
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Centrality dependence is implicit in Eq. J7J and the mod- 
ified cross section per average collision da^ A ^ h " / (N' t 



coll \ 
dAul 



can be calculated from Eqs. J2J and © via the substi- 
tution Eq. ©. The necessary scaling with impact pa- 
rameter, b, of the enhanced power corrections can be ob- 
tained through the nuclear thickness function, A 1 / 3 — » 
A 1 /3T j4 (6)/T A (6 roin . bias ), where T A (b) = p A (z.b)dz. 

The top panels of Fig. [21 show the rapidity and trans- 
verse momentum dependence of R^ Au (b) and R^ A ^(b) 
for minimum bias collisions. The amplification of the 
suppression effect at forward y\ comes from the build-up 
of coherence and the decrease of the Mandelstam vari- 
able (— t). At high pt\iPti the attenuation is found to 
disappears in accord with the QCD factorization theo- 
rems jy, [lj- BRAHMS data [j] on h~ production at 
forward rapidity is shown for comparison. Our result 
should be contrasted with a recent numerical study of 
non-linear parton evolution that found a factor of 2 sup- 
pression of the hadron production in p+A collisions at 
Pt = 10 2 — 10 4 GeV but negligible rapidity dependence 
of Rp A (b) at small p T < 2 GeV The bottom pan- 

els of Fig. [21 show the growth of the nuclear attenuation 
effect with centrality. 

In the double collincar approximation, Eq. J3J), the 
leading hadrons are produced strictly back-to-back. The 
acoplanarity, Atp ^ tt |l9j| . arises from soft gluon resum- 
mation [2(j , next-to- leading order corrections and in the 
presence of nuclear matter - transverse momentum diffu- 
sion [l4|. We consider dihadron correlations associated 
with hard scattering partonic 2 — > 2 processes that have 
the bulk many-body collision background subtracted and 
can be approximated by near-side and away-side Gaus- 
sians: 



C 2 (Acp) 



1 



" Jv dijct 



Norm dAip 



27T(7Nc 



+ 



iFar 



27T(7p a 



(8) 



In Eq. © the near-side width ONear of C2(Aip) is deter- 
mined by jet fragmentation 21] . If the strength of the 
away-side correlation function in elementary N+N colli- 
sions is normalized to unity, dynamical quark and gluon 
shadowing in p+A reactions will be manifest in the at- 



tenuation of the 
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L pA '-(b). The away-side 
width (7F al . is related 0, 0] to the accumulated dijet 
transverse momentum squared in a plane perpendicular 
to the collision axis: 
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In Eq. yi are the interacting parton rapidities and 
/j, 2 L e ff/X = 0.72 GeV 2 characterizes the elastic scat- 
tering power of cold nuclear matter in minimum bias 
d+Au collisions 0, 0] and varies with centrality ~ 
T A (b)/T A (b min . hias ). The vacuum broadening (k T ) mc is 
taken from the data [2ll |22| 

In triggering on a large transverse momentum hadron, 
the nuclear modification from enhanced power correc- 
tions and transverse momentum diffusion is reflected in 
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FIG. 3: Left: centrality dependence of CzfA^?) at moderate 
PTi = 4 GeV, pt 2 = 2 GeV and rapidities yi = 0, 2 and 
y- 2 = 0. Central d+Au and p+p data from STAR 22]. Right: 
C2(A</?) at small transverse momenta pr x = 1.5 GeV, pr 2 
1 GeV and rapidities yi = 0, 4 and 3/2 = 0. 



the away-side correlation function. For large vacuum 
broadening, (k^) vac = 2.5 — 3.5 GeV 2 , elastic scatter- 
ing may lead to only a moderate additional growth of 
(TFar- The left panels of Fig.[3]show that for pr 1 = 4 GeV, 
Pt 2 = 2 GeV the dominant effect in C2 (A</j) is a small in- 
crease of the broadening with centrality, compatible with 
the PHENIX [2l[ and STAR ^ measurements. Even at 
forward rapidity, such as y% = 2, the effect of power cor- 
rections in this transverse momentum range is not very 
significant. At small — 1.5 GeV, px 2 = 1 GeV, shown 
in the right hand side of Fig. |3J the apparent width of 
the away-side ^(At/j) is larger. In going from midrapid- 
ity, yi = j/2 = 0, to forward rapidity, y x = 4, y 2 = 0, 
we find a significant reduction by a factor of 3 - 4 in the 
strength of dihadron correlations from the nuclear en- 



hanced power corrections. We emphasize again that this 
result depends sensitively on centrality and the choice of 
transverse momenta and disappears at high px- 

The dynamical cross section attenuation calculated 
here does not contradict, instead, complements the ef- 
fects from a possible modification of the nuclear parton 
distribution functions (nPDFs), known as leading twist 
shadowing 23] . In our formalism, one can include both 
effects by replacing the PDFs in Eqs. @ and © with 
the corresponding nPDFs and applying Eq. ©. The 
weaker scale dependence of nPDFs will thus slow down 
the disappearance of the nuclear suppression at high px 
in Fig. [21 A combined global QCD analysis would, how- 
ever, be required and is beyond the scope of this Letter. 

In conclusion, we presented a systematic approach to 
the calculation of coherent QCD multiple scattering and 
resummed the nuclear enhanced power correction to the 
single and double inclusive hadron production cross sec- 
tions in p+A reactions. At low pr we find a sizable sup- 
pression, which grows with rapidity and centrality. At 
high pt the nuclear modification disappears in accord 
with the QCD factorization theorems. We demonstrated 
that both particle spectra and dihadron correlations are 
sensitive measures of such dynamical attenuation of the 
parton interaction rates. 

Our approach, with its intuitive and transparent re- 
sults, can be easily applied to study the nuclear modi- 
fication of other physical observables in p+A reactions 
and its predictions can be readily tested against RHIC 
data. The systematic incorporation of coherent power 
corrections provides a tool to address the most interest- 
ing transition region between "hard" and "soft" physics 
in hadron- nucleus collisions. It allows for the extension of 
perturbative QCD calculations to relatively small trans- 
verse momenta and for bridging the gap between the par- 
ton model picture and the possible onset of gluon satura- 
tion [2j, |24j at very large collider energies and very small 
values of nuclear x. 
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